Prom caging to Rouse dynamics in polymer melts with intramolecular barriers: 
a critical test of the Mode Coupling Theory 
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By means of computer simulations and solution of the equations of the Mode Coupling Theory 
(MCT), we investigate the role of the intramolecular barriers on several dynamic aspects of non- 
entangled polymers. The investigated dynamic range extends from the caging regime characteristic 
of glass-formers to the relaxation of the chain Rouse modes. We review our recent work on this 
question, provide new results and critically discuss the limitations of the theory. Solutions of the 
MCT for the structural relaxation reproduce qualitative trends of simulations for weak and moderate 
barriers. However a progressive discrepancy is revealed as the limit of stiff chains is approached. This 
disagreement does not seem related with dynamic heterogeneities, which indeed are not enhanced 
by increasing barrier strength. It is not connected either with the breakdown of the convolution 
approximation for three-point static correlations, which retains its validity for stiff chains. These 
findings suggest the need of an improvement of the MCT equations for polymer melts. Concerning 
the relaxation of the chain degrees of freedom, MCT provides a microscopic basis for time scales from 
chain reorientation down to the caging regime. It rationalizes, from first principles, the observed 
devations from the Rouse model on increasing the barrier strength. These include anomalous scaling 
of relaxation times, long-time plateaux, and non-monotonous wavelength dependence of the mode 
correlators. 

PACS numbers: 64.70.pj, 64.70.qj, 61.20.Ja, 83.10.Mj, 83.10.Rs, 83.80.Sg 
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I. INTRODUCTION 

The different dynamic processes present in amorphous 
polymers cover a extremely broad range of characteristic 
time scales, spanning from about 100 femtoseconds up 
to years. There are two main reasons for this. First, 
polymers are usually good glass-formers, which inher- 
ently exhibit a dramatic increase of the viscosity and 
structural (a-) relaxation times on approaching the glass 
transition temperature Tg. As in non-polymeric glass- 
formers, localized dynamic processes are also present be- 
low Tg pj]. Second, their macromolecular character in- 
troduces relaxation processes related to the dynamics of 
the internal chain degrees of freedom. In the case of low- 
molecular weight, nonentangled, polymer chains a sub- 
linear increase (Rouse-like) arises in the mean squared 
displacement prior to the linear diffusive regime. In the 
case of high-molecular weight, strongly entangled, chains 
further sublinear regimes are found between the Rouse 
and linear regimes, which are usually interpreted in terms 
of reptation dynamics [2l-fi] . Such processes are inherent 
to chain connectivity, and extend over more time decades 
on increasing chain length. This broad time window for 
chain dynamics is observed even for temperatures far 
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above Tg, when the structural relaxation extends over 
just a few picoseconds. 

Another particular ingredient of polymer systems is 
that, apart from fast librations or methyl group rotations 
0, every motion involves jumps over carbon-carbon ro- 
tational barriers and/or chain conformational changes. 
The corresponding map of relaxation processes is largely 
influenced by the barrier strength. Intramolecular barri- 
ers play a decisive role injC.g., crystallization [1, ad- 
sorption onto surfaces [1, Q , viscoelastic properties [To| , 
or phase behavior of block copolymers [ll| . Models for 
semiflexible and stiff polymers are of great interest in 
biophysics, since they can be applied to many impor- 
tant bi opolym ers as DNA, rodlike viruses, or actin fila- 
ments |12h14| . Thus, an understandig of the role of the 
intramolecular barriers on structural and dynamic prop- 
erties of polymer systems is of practical as well as of 
fundamental interest in many fields of research. 

A possible theroretical approach to this problem is 
provided by the Mode Coupling Theory (MCT) 0. 
MCT introduces a closed set of coupled Mori-Zwanzig 
equations for the time dependence of density correlators. 
Static correlations enter the memory kernel as external 
input. Since the former can be related to the interaction 
potential through liquid state theories, MCT constitutes 
a first-principle theory for slow dynamics in complex sys- 
tems. MCT has been developed over the last years to 
include systems with intramolecular structure (see e.g.. 
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[iGl - fTsj V This includes the approach of Chong and co- 
workers for simple polymer melts [l^ [l^ , based on the 
polymer reference interaction site model (PRISM) [2l| 
for the static correlations. This approach was applied 
to the specific case of fully-flexible chains t20i], i.e., 
without intramolecular barriers. A major success was 
the derivation, from first-principles, of the scaling laws 
predicted by the phenomenological Rouse model [2| for 
chain dynamics in nonentangled polymer melts. Like- 
wise, it provided a unified microscopic description of both 
chain dynamics and the structural relaxation associated 
to the glass transition [l^, [2^ . 

Some of us have recently performed a systematic com- 
putational investigation of the role of intramolecular bar- 
riers on the glass transition in polymer systems [13, . 
Starting from fully-flexible bead-spring chains, we intro- 
duced stiffness by implementing intramolecular barriers 
with tunable bending and torsion terms. In Ref. |23| we 
discussed the glass transition within the framework of the 
MCT for polymer melts, comparing simulations with nu- 
merical solutions of the MCT equations, in the long-time 
limit, for a broad range of barrier strength. This was 
possible since the quality of the PRISM approximations 
observed for fully-flexible chains |2^ was not affected at 
all by the introduction of internal barriers in all the in- 
vestigated range [23| . Numerical solutions reproduced 
trends in the nonergodicity parameters and MCT critical 
temperatures for weak and moderate barriers. However, 
strong discrepancies were observed on approaching the 
limit of stiff chains p^ . 

In this article we briefly summarize the main points of 
Refs. [l^ll^ and present extensive new results. Thus, we 
solve the time-dependent MCT equations for density cor- 
relators, and compare simulation and theoretical trends 
in a-relaxation times. We critically discuss the limita- 
tions of the theory by analyzing the accuracy of the as- 
sumed approximations. We find that dynamic hetero- 
genities, static three-point correlations, and chain pack- 
ing effects not accounted by MCT, do not play a major 
role on increasing the barrier strength. Indeed their ef- 
fects seem to be weaker that in the case of fully-flexible 
chains. The reason for the observed discrepancies be- 
tween simulation and theory for very stiff chains remains 
to be understood. 

We also present here a systematic investigation on the 
effect of intramolecular barriers on the internal chain dy- 
namics of nonentangled polymers. We analyze correla- 
tors for the chain normal modes (Rouse modes) and for 
bond reorientation. The simulations reveal strong devia- 
tions from the Rouse model on increasing chain stiffness. 
These include anomalous scaling of relaxation times [1^ , 
long-time plateaux, and nonmonotonous wavelength de- 
pendence of the mode correlators. We show that these 
anomalous dynamic features are reproduced by the cor- 
responding MCT equations for the Rouse modes. This 
generalizes the analysis of Ref. [1^, which was limited 
to fully-flexible chains, to polymers with intramolecular 
barriers of arbitrary strength. Thus, beyond usual phe- 



nomenological models for chain dynamics, MCT provides 
a unified microscopic picture down to time scales around 
and before the a-process [1^ . 

The article is organized as follows. In Section II we de- 
scribe the model and give simulation details. In Section 
III we compare simulation results with MCT solutions 
for several dynamic correlators probing structural relax- 
ation and chain dynamics. In Section IV we discuss the 
possible origin of the observed deviations from MCT pre- 
dictions. Conclusions are given in Section V. 

II. MODEL AND SIMULATION DETAILS 

We have performed molecular dynamics (MD) simu- 
lations of a bead-spring model with tunable intramolec- 
ular barriers. All chains consist of N^^ — 10 identical 
monomers of mass m = 1. Non-bonded interactions be- 
tween monomers are given by a corrected soft-sphere po- 
tential 

V{r)^4e[{<7/ry'-a + C2{r/an (1) 

where e — 1 and a = 1. The potential V{r) is set to 
zero for r > ca, with c = 1.15. The values Cq = 7c~^^ 
and C2 = 6c~^'^ guarantee continuity of potential and 
forces at the cutoff distance r ^ ca. The potential V{r) 
is purely repulsive. It does not show local minima within 
the interaction range r < ca. Thus, it drives dynamic 
arrest only through packing effects. Chain connectivity 
is introduced by means of a finitely-extensible nonlin- 
ear elastic (FENE) potential ,28,] between consecutive 
monomers: 

^FENE(r) = -eKpRl ln[l - {Roa)-\% (2) 

where Kp — 15 and Rq — 1.5. The superposition of 
potentials ^ and © yields an effective bond potential 
for consecutive monomers with a sharp minimum at r « 

0. 985, which makes bond crossing impossible. 
Intramolecular barriers are implemented by means of 

the combined bending and torsional potentials proposed 
by Bulacu and van der Giessen in Refs. [2^, Is^l- The 
bending potential Vb acts on three consecutive monomers 
along the chain and is defined as 

VB{e,) = (ei^B/2)(cos^,; - COS^o)', (3) 

where 6i is the bending angle between consecutive 
monomers i — 1, i and i + 1 (with 2 < i < N^i — 1). 
We use ^0 — 109.5° for the equilibrium bending angle. 
The torsional potential Vt constrains the dihedral angle 
ipi^i+i. The latter is defined for the consecutive monomers 
i-l, i, i + l and i + 2 (with 2 < « < iVm - 2), as the 
angle between the two planes defined by the sets (i — 1, 

1, i + \) and (i, i + 1, i + 2). The form of the torsional 
potential is 

Ft(0„0.+i,0.,.+i) = 

3 

eKt sin^ 9i sir? 9i+i ^ a„ cos" (f>i,t+i . (4) 

n=0 
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The values of the coefficients a„ are on = 3.00, oi = 
-5.90, 02 = 2.06, and 03 = 10.95 [11,130]. The torsional 
potential depends both on the dihedral angle and 
on the bending angles 6i and 6*2+1. As noted in Refs. jH, 
[13 , the functional form Q avoids numerical instabilities 
arising when two consecutive bonds align, without the 
need of imposing rigid constraints on the bending angles. 

In the following, temperature T, time t, distance, 
wave vector q, and monomer density p are given respec- 
tively in units of e/ke (with the Boltzmann con- 
stant), CT(m/e)^/^, a, a^^, and a^^. We investigate, at 
fixed monomer density p — 1.0, the temperature depen- 
dence of the dynamics for different values of the bending 
and torsion strength, {Kb,Kt) = (0,0), (4,0.1), (8,0.2), 
(15,0.5), (25,1), (25,4), and (35,4), covering a broad dy- 
namic range from the caging characteristic time to the 
relaxation time of the slowest Rouse mode. We inves- 
tigate typically 8-10 different temperatures for each set 
of values (Kb,Kt)- Additional numerical details can be 
found in Refs. [11, Hlj. 

The investigated range of barrier strength corresponds 
to a strong variation of the chain stiffness. This can be 
quantified by the average end-to-end radii, Ree, of the 
chains. Thus, for the representative values {Kb,Kt) = 
(0,0), (8,0.2), (25,1) and (35,4), which cover the range 
from fully-flexible chains to the stiffest investigated 
chains, we find i?oo — 3.6, 4.7, 5.5, and 6.5 at the re- 
spective lowest investigated temperature. 

III. RESULTS: SIMULATIONS vs. THEORY 

Before addressing the dynamic aspects of the present 
system, we want to stress that the investigated state 
points correspond to isotropic phases. We do not observe 
signatures of global orientational order induced by chain 
stiffness for the investigated state points. Thus, by mea- 
suring the quantity ^2(6) = (3(cos^ 6) — l)/2, where Q 
is the angle between the end-to-end vectors of two chains, 
and averaging it over all pairs of distinct chains, we ob- 
tain in all cases values |P2(0)| < 3 x lO^'^. This is illus- 
trated in Fig.[I] which shows the time evolution of P2(0) 
along a typical simulation window, both for fully-fiexible 
chains, {Kb,Kt) = (0,0), and for the stiffest investigated 
chains, IKb,Kt) = (35,4). 

Local orientational order is also negligible. This is evi- 
denced by computing a similar correlator J-2(0; ?'cm)- In 
this case the average is performed only over pairs of dis- 
tinct chains for which the distance between their respec- 
tive centers-of-niass is less than rem- Fig. [T] displays, for 
the former cases of fully-flexible and stiff chains, data of 
-P2(0; rem) for several values of rem. Negligible values of 
^2(0; fcm) are obtained for rem > 2.0. Thus, the time av- 
erage over the simulation time window, tgim, provides val- 
ues |(P2(0;?'cm > 2.0))t,.„| < 0.02. By comparing both 
panels we conclude that chain stiffness does not induce a 
significant increase, if any, of local orientational order in 
the investigated systems. Weak local orientational order 
|(P2(0; ?'cm))tai,„ I ^ 0.1 is observed only for very small in- 
terchain distances (see data for P2(0; ''cm = 1-4)). Again, 
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FIG. 1: Time evolution of the global and local orientational 
parameter (see text) at two selected temperatures, for fully- 
flexible chains (a) and stiff chains (b) with {Kb,Kt) ~ (35, 4). 



the introduction of chain stiffness does not induce clear 
changes in the orientational order at this length scale. 

A. Structural relELxation 

Now we characterize dynamic features associated 
to the caging regime and the structural a-relaxation. 
Fig. [2^ shows the self-density correlator f^{q,t) at fixed 
T = 1.5 and for several values of the barrier strength. 
The former is defined as P{q,t) — -/V~^(X)jLi exp[iq • 
{rj{t) — rj{0))]). The sum is done over the coordi- 
nates Tj of all the N monomers in the system. In 
all the cases the correlator is evaluated at the max- 
imum, (/max ~ 7 [23], of the static structure factor 

S{q) = N-\j:f,k^i^M^^ ■ (r,(0) - rfc(O))]). We ob- 
serve that increasing the strenght of the internal bar- 
riers at fixed p and T leads to slower dynamics. In 
the fully fiexible case f^{q,t) decays to zero in a single 
step. On increasing the strength of the internal barriers 
f^{q,t) exhibits the standard behavior in the proximity 
of a glass transition. After the initial transient regime, 
P{q,t) shows a first decay to a plateau, which is associ- 
ated to the caging regime, i.e., the temporary trapping 
of each particle by its neighbors. At long times, a second 
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FIG. 2: Panel (a): self-density correlator /''(gmax, i), for 
several values of the barrier strength at fixed temperature 
T = 1.5. Panel (b): temperature dependence of the relax- 
ation times To. 2 of /"(gmaxjt), for the former values of the 
barrier strength. 



decay is observed from the plateau to zero. This corre- 
sponds to the structural a-relaxation. Similar trends are 
displayed by the density-density correlator (not shown), 
defined as fiq,t) = (p(q, t)p(-q, 0))/(p(q, 0)p(-q, 0)), 
with p{ci,t) = Ej^iexp[iq- rj(t)]. 

Let us define to.2 as the time for which /*'(gmax, ''0.2) — 
0.2, and r?- as that obtained from fitting the a- 



decay to a Kohlrausch- Williams- Watts (KWW) function, 
Agexp[-{t/T^)l^] (with Aq,l3 < 1). Both times cor- 
respond to a significant decay from the plateau, and 
therefore can be used as operational definitions of the a- 
relaxation time Tq,. Fig. ^jp shows T0.2 as a funcion of T, 
for different values of the bending and torsional constants 
(results for are analogous). As observed in the anal- 
ysis of the self-correlators, increasing the chain stiffness 
slows down the dynamics. At fixed temperature, the re- 
laxation time for the stifFest investigated chains increases 
by several decades with respect to the fully-flexible case. 

The dynamic trends displayed in Fig. [5] demonstrate 
that intramolecular barriers constitute an additional 



mechanism for dynamic arrest, coexisting with the gen- 
eral packing effects induced by density and temperature. 
Now we discuss this scenario within the framework of the 
(ideal) MCT. We briefly summarize the basic concepts 
and predictions of the theory. Extensive reviews can be 
found, e.g., in Refs. [Tsl lsil - lssj . On approaching a glass 
transition from the ergodic phase, density fluctuations 
decay in a slower fashion, remaining frozen in amorphous 
configurations when the glass transition occurs. MCT de- 
scribes this phenomenon as a feedback mechanism driven 
by the slow density fluctuations. By starting from the 
fundamental Liouville equation of motion and using the 
Mori-Zwanzig projection operator formalism, an integro- 
differential equation is obtained for the density-density 
correlator: 



f{q,t) 



mS{q) 



mS{q) ' 

dt'm{q,t -t')f{q,t') = 0. 



(5) 



The memory kernel m{q,t - t') cx (i?^(0)i?^(t - t')), is 
expressed in terms of the associated fluctuating forces 
i?q In order to provide a closed solvable form of 

Eq. ([5]), MCT introduces several approximations for the 
memory kernel. These approximations are: 

i) The fluctuating force can be splitted in two terms: the 
regular (fast) contribution, linear in density fluctuations, 
and a second term which can be expressed as a a 
linear combination of 'mode pairs', pkPq-k- The latter 
provides the slow contribution relevant for the structural 
relaxation. Thus, the first approximation consists of 
neglecting the fast contribution. 

ii) Convolution approximation: three-point static 
correlations are approximated as products of static 
structure factors. 



(p_q(0)pu(0)pq_k(0)) « NS{q)S{k)S{\ci - k|). (6) 

iii) Kawasaki approximation: dynamic four-point corre- 
lations are factorized in terms of products of dynamic 
two-point correlations (see e.g., Ref. |33] for details). 
Nowadays there is plenty of evidence that this approx- 
imation worsens on decreasing temperature, specially 
around the time scale of the a-relaxation. The break- 
down of the former approximation is usually assigned to 
the emergence of strong dynamic heterogeneities in the 
proximity of the glass transition [36l - l40j . 

After applying the former approximations, the memory 
kernel m{q,t) becomes bilinear in f{q,t): 



■m{q,t) 



(2^ 



V(q,q-k)/(fc,t)/(|q-k|,i). (7) 
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The vertex V(q, q — k) is given by: 

V(q,q-k)- 2^5(q)5(fc)5(|q-k|)x 

[q•kc(fc) + q•(q-k)c(|q-k|)]^ (8) 

where c{q) is the direct correlation function [ilj. Eq. ^ 
constitutes a closed set of coupled equations which can 
be solved self-consistently, provided S{q) and c{q) are 
known. The latter are external inputs in the MCT equa- 
tions. Since static correlators contained in the vertex 
vary with the control parameters (e.g., density, tempera- 
ture, or barrier strength), the MCT equations ([7]) estab- 
lish a direct connection between statics and dynamics. 
Moreover, the former static correlators can be related to 
the interaction potential through closure relations from 
liquid state theories [41]. With this, MCT provides a 
first-principle approach for the slow relaxation of density 
correlators. 

Recently, Chong and co-workers have derived MCT 
equations for simple models of polymer melts [lol [20| . By 
exploiting the polymer reference interaction site model 
(PRISM) ^2lj, the MCT equations are considerable sim- 
plified. This is achieved by replacing site-specific inter- 
molecular surroundings of a monomer by an averaged 
one (equivalent site approximation), whereas the full in- 
tramolecular dependence is retained in the MCT equa- 
tions 1&, 20]. The so-obtained scalar MCT equations of 
motion, memory kernel and vertex for polymer chains are 
formally irfeniica/ to Eqs. (|5l7l8p . The polymer character 
of the system only enters implicitly through the PRISM 
relation [2l| pc{q) = \/uj{q) — l/Sjq), which differs from 
the Ornstein-Zernike equation [41|, pc{q) = 1 — S^^{q), 
for monoatomic systems. The quantity uj{q) is the chain 
form factor, defined as 

^(9) - E E (e^p[*q • (^a(o) - rim]) , (9) 

" 1=1 a, 6=1 

where are the coordinates of the ath monomer in the 
Ith chain. A^c is the total number of chains. The use of 
the former MCT equations is a priori justified for poly- 
mers of variable stiffness. Indeed, it has been shown that 
the PRISM approximations retain their validity not only 
in the fully-flexible limit [l^li but also when strong in- 
tramolecular barriers are present [23l |. 

For the case of the self-density correlators P{q,t), the 
MCT equations are different from the monoatomic case. 
The former are obtained by summation of the diagonal 
terms of the self site-site density correlators. The latter 
are given by 

1=1 

with indices defined as in Eq. ©. The correlators 
F^fj{q,t) are determined by sol ving the corresponding 
MCT matrix equations (see Ref. j20|). 



Ideal MCT predicts a sharp transition from an er- 
godic liquid to an arrested state (glass), at a given value 
of the relevant control parameter (temperature in the 
present case). At the transition (or 'critical') tempera- 
ture T — Tc, the non-ergodicity parameter, defined as 
fq — lim4_).oo /(<Z,i), jumps from zero to a nonzero value 
fq. The latter is called the critical non-ergodicity param- 
eter. By taking the limit t — >■ oo in the MCT equations, 
one finds the relation 

Eq. pT|) always has the trivial solution {fq} = 0. Glassy 
states take place when solutions fq > also exist. The 
temperature at which the jump from zero to nonzero so- 
lutions occurs defines Tc. The corresponding solutions 
define the critical non-ergodicity parameters. 

The separation parameter, ex = [T — Tc)/Tc measures 
the distance to the critical temperature. We are inter- 
ested in the behavior of f{q,t) in the ergodic fiuid, i.e., 
for > 0. For small values of ct, MCT predicts several 
asymptotic laws for dynamic observables [1^, which are 
characterized by several dynamic exponents. Thus, the 
decay from the plateau follows a von Schweidler expan- 
sion f{q,t) = f^ - hqt^ + h^qh"^^ + 0{t^''), and diffusivi- 
ties and a-relaxation times obey D^^ , Ta ^ (T — Tc)^^ . 
The exponents of these asymptotic laws are related to 
the so-called exponent parameter A, which is the only 
independent one. The MCT expression for A is deter- 
mined by the static correlators S{q) and c{q) evaluated 
at T = Tc, and b y th e critical non-ergodicity parameters 
(see, e.g., ^M). 

We solved Eqs. ([5]) and dTTI) by combining simula- 
tion results of uj{q) with the PRISM equation pc{q) — 
l/ui(q) — 1/S{q) and the Percus-Yevick closure relation 
[4l| . Details of the numerical procedure for solving (fTTj) 
can be found in Ref. [23i]. Numerical integration of the 
density correlators was performed following the method 
of Ref. [i^l. It often happens in the analysis of exper- 
iments or simulations that numerical solutions of the 
MCT equations are not available. In such cases, a phe- 
nomenological analysis can be performed, and the values 
of Tc and the former dynamic exponents can be obtained 
as fit parameters from the experimental or simulation 
data. Consistency of the analysis requires that the expo- 
nents, which are obtained from independent fits to differ- 
ent scaling laws, are related to the same A-parameter, as 
predicted by the theory. This consistency test was done 
in the analysis of our simulation data (see Refs. [l^, 
for a detailed explanation), providing different values of 
Tc and A for each barrier strength. These values obtained 
from simulations can be compared with the values pro- 
vided by solution of the MCT equations. 

This comparison is shown in Fig. [31 Superscripts 'MD' 
and 'MCT' are used respectively for simulation and the- 
oretical values. The data are represented as a function 
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FIG. 3: Critical temperature Tc (a) and A-exponent (b) as 
a function of the end-to-end radius R^^ (measured at T^^). 
The theoretical values T^'~^'^ and A^*-'^ are compared with the 
simulation values and (see text). The respective 

values of the bending and torsional strength {Kb,Kt) are 
indicated on top of each symbol for in (a) and for A^"^ 

in (b). 



FIG. 4: For fully-flexible and stiffest investigated chains, q- 
dependence of the KWW times for self-density correlators 
p{q,t). Data are rescaled by the respective KWW times 
at gmax- Data in panels (a) and (b) correspond respectively 
to simulation results and MCT solutions at selected temper- 
atures (see legend). 



of the end-to-end radius, which quantifies chain stiff- 
ness. A clear correlation between the barrier strength 
and the values of and A^'^ is unambiguously demon- 
strated. The interplay between monomer packing effects 
and intramolecular barriers j2S] induces a progressive in- 
crease of at fixed density. We note that from the 
fully-flexible limit {K^,Kt) = (0,0) to barriers with 
{Kb,Kt) = (15,0.5), the data sets for T^^ and T^^'^ 
roughly display the same slope. As usual, there is a shift 
factor between simulation and theoretical temperatures 
(here yMD/j^MCT _ i 25), which may originate from the 
mean-field character of MCT [l^ . The range of barrier 
strength for which J'^'^CT ^^^^^^ are roughly parallel 
is significant. Indeed, for {Kb,Kt) = (8,0.2) the end- 
to-end radius R^^ is a 30% larger than for fully-flexible 
chains. However, a strong discrepancy between simula- 
tion and theory becomes evident on increasing the barrier 
strength from {Kb,Kt) = (15,0.5). While beyond this 



point T^'-^'^ seems to approach an asymptotic limit, 
increases up to 1.23 for the stiffest investigated chains. 

Similar trends are observed for the A-exponent. Simu- 
lation values increase from A^^ = 0.76 for fully-flexible 
chains to A'^^ = 0.86 for the stiffest investigated chains. 
The first ones are typical of simple glass-formers as the 
archetype hard-sphere fluid (A = 0.74 J^]), where dy- 
namic arrest is driven by packing effects. The largest ones 
are similar to those observed in realistic models of poly- 
mer melts which incorporate the full chemical structure 
of the chains [43 - |47j . On the contrary, the theoretical 
exponent exhibits a very weak variation, 0.71 < A'^^"'" < 
0.72, over the investigated range of barrier strength. 

These discrepancies in the case of strong intramolec- 
ular barriers are also reflected in the g-dependence of 
density correlators computed from simulations and from 
solution of the MCT equations. In both cases we flt- 
ted the corresponding a-decay to a KWW function (see 
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above). Fig. |4] compares the g-dependence, at fixed T, 
of the KWW time for the self-correlators f^{q,t), as 
obtained from simulations and from theory. Results are 
presented for the fully-flexible case and for the stifFest 
investigated chains. 

Before discussing such results, some points must be 
clarified. As mentioned above, the mean-field character 
of MCT usually yields a temperature shift between simu- 
lation and theory (see Fig.|3^). Moreover, MCT times are 
affected by an undetermined constant factor [l^. Thus, 
a proper comparison between theory and simulation for 
time-dependent correlators can be done by rescaling t by 
some characteristic relaxation time, and using a common 
separation parameter er [lOl- This is the case for the 
data of Fig. U) Thus, each data set is rescaled by the re- 
spective KWW time t^„^^ corresponding to /'^(<7max, i)- 
Temperatures of both panels correspond to et ~ 0.04 
and 0.08 for respectively fully-fiexible and stiff chains. In 
the rest of the article we will present several comparisons 
between simulation data and MCT solutions. It will be 
understood that the times and temperatures of the com- 
pared data obey the former criteria. 

A clear disagreement between simulation and theoret- 
ical trends becomes evident in Fig. 31 The two sets of 
KWW times obtained from simulation show a rather dif- 
ferent g-dependence, which is more pronounced for the 
stiff chains. On the contrary, after rescaling by Tq^^^, 
the theoretical sets become esentially identical. Fig [S] 
shows a similar comparison between simulation and the- 
ory for the rescaled KWW times of the density-density 
correlators, f{q,t). For q > Qmax the theory reproduces 
qualitatively the shape of the relaxation times, which are 
modulated by the respective static structure factor S{q) 
(not shown). However, MCT fails at reproducing the 
broad peak at intermediate ~ 4 which is present in 
the simulation data. This failure was already noted for 
fully-flexible chains in Ref. (2^, and is confirmed here 
for the general case with intramolecular barriers. Ap- 
parently (note the error bars), the peak does not shift 
significantly and decreases its intensity as chains become 
stiffer, leading to a shoulder. In previous works [i^ on 
similar fully-flexible bead-spring chains of — 10, the 
value of qc has been identified with 27r/i?g, where i?g is 
the chain radius of gyration. Data of Fig. [S^ do not seem 
compatible with this assignment. Appart from results for 
fully-flexible and stiffest investigated chains of = 10, 
we include data for {Kb,Kt) = (15,0.5) of additional 
simulations with N^n = 21. With this, the data sets of 
Fig. [5^ cover a significant variation in i?g. Namely, for 
{Kb,Kt) = (0,0), (35,4), and (15,0.5) we respectively 
find 27r/i?g = 4.2, 2.8, and 2.0. Thus, the observation 
qc ~ 27r/i?g for fully-fiexible chains is apparently fortu- 
itous. The associated length scale 27r/gc ~ 1.6<J rather 
seems to be a characteristic feature which does not de- 
pend significantly on the barrier strength. We will come 
back to these points in Section IVc. 

B. Chain dynamics 
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FIG. 5: Main panels: for fully-flexible and stifl^est investigated 
chains, g-dependence of the KWW times for density-density 
correlators f{q,t). Data are rescaled by the respective KWW 
times at Qmax. Data in (a) and (b) correspond respectively 
to simulation results and MCT solutions at selected temper- 
atures (see legend). Typical error bars are included in (a). 
Dashed lines are guides for the eyes. Inset of panel (a): Static 
structure factors for the chains centers-of-mass, Scmiq)- 



In this subsection we compare simulation and theo- 
retical results for the dynamics of the Rouse modes. 
First we briefly summarize the assumptions and main 
predictions of the Rouse model. The starting point 
is a tagged gaussian chain of A^„i monomers connected 
by harmonic springs of constant S/cbT/^^, with b the 
bond length. The effective interaction experienced by 
the monomers is given by a friction coefficient C ^^id 
a set of stochastic forces ij. Excluded volume inter- 
actions are neglected. The chain motion is mapped 
onto a set of A^m normal modes (Rouse modes) la- 
belled by p = 0, 1, 2, A^ni — 1, of wavelength N/p, 
and defined as 0, |4§| Xp(i) = J2f=i Pjp^ji't)^ '^i^h 
P3P = V(2-'5po)/A^mCOs[(j - l/2)p7r/A,J. The chain 



8 



0.05 ■ 



(a) 



(Kg,K^) = (0,0) 



• T=0.50 
■ T=1.50 



o 

a 



10 



o 



MD 



. (K3,K,) = 


(0,0) 


T = 


0.55 


. (K3,K,) = 


(8,0.2) 


T = 


0.70 


. (K3,Kp = 


(25,1) 


T = 


1.06 



-0.05 ■ 
0.6 ; 

'k 0.4 ; 

I 

0.2 ; 



I 

lo" 



(b) 



(K3,K^) = (35,4) 



— I 

T= 1.22 
T= 1.48 
T= 1.90 



10' 







20 



40 



60 



9(p-l) + q 



80 



(c) 



<U 



(K^,K.p = (0,0) T = 
(Kg,K.p = (0,0) T=l. 
(Kg.ICj) = (8,0.2) T = 
(Kg,K,) = (8,0.2) T = 



50 
50 
0.64 
1.50 



10 r 



o 







A (Kg,IC,) = (25,1) 


T = 


0.96 


10-^ 




^ (K3,K.j) = (25,l) 


T = 


1.75 






. (Kg,K^) = (35,4) 


T = 


1.22 






□ (K3,IC,) = (35,4) 


T = 


1.90 



H □ 



FIG. 6: Static intrachain correlations computed from simula- 
tions. For each value of {Kb, Kt) (see legends) results include 
data at the highest and lowest investigated T. Panels (a) and 
(b): off-diagonal terms of *p<j(0) (see text) versus the com- 
bined variable 9{p ~ 1) + q. Data in (a) and (b) respectively 
correspond to fully-flexible chains and to the stiffest investi- 
gated chains. Panel (c): diagonal terms Cpp{0) (see text). 
Each data set corresponds to a fixed value of {Kb,Kt) and 
T (see legend). For clarity, each set has been rescaled by a 
factor a in order to get a(7ii(0) = 1 in all cases. Dashed lines 
indicate approximate power-law behavior ~ p^^- From top 
to bottom, X = 2.2, 2.7, 3.1, and 3.8. 



center-of-mass coincides vi^ith J^o{t)/^/N^. The mode 
correlators are defined as Cpq{t) = [(Xp(0) • Xg(t)) — 
<5o,px9(Xp(0) • Xg{0))]/3N^. For p,q > we define 
tfie matrix C'pg(O) = Cpg(O). In the Rouse model 
the stochastic forces are fully spatial and time uncor- 
related, i.e., {fj{t) ■ (kit')) = 6CkBT6jkS{t - t'). These 
properties of the random forces lead to orthogonal- 
ity and exponentiality of the Rouse modes [2|. Thus, 
the mode correlators obey Cpq{t) = C'pg(O) exp[— t/xp], 
with Cpq{Q) ^ ^pg(6V24iV^)sin"2[p7r/2iVm] and Tp = 
(C&^/12A;BT)sin-2[p7r/2iVm]. Accordingly, for p < iV,„ 
the quantities C'pj,(0) and Tp scale as ~ 

In the following we show how the former scaling prop- 
erties are strongly altered by the introduction of in- 
tramolecular barriers. This is demonstrated in Fig. [6] 
for the case of intrachain static correlations. We 
show the off-diagonal terms of 'i!pq{Q) — (Xp(0) • 
Xg(0)/(Xp(0)Xq(0))) (the diagonal terms are trivially 
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FIG. 7: Simulation results (filled symbols) and MCT solutions 
(empty symbols) , for the p-dependence of the relaxation times 
Tp of the mode correlators. Each data set corresponds to a 
fixed value of {Kb, Kt) and T (see legends). For clarity, each 
set is rescaled by its respective t\. Dashed and dotted lines 
indicate approximate power-law behavior ~ ■ From top 
to bottom, simulations (dashed): x = 2.0, 2.4, 2.8; MCT 
(dotted): x = 2.2, 2.7, 3.1. 



^'pp(O) = 1). Data for fully-flexible chains exhibit small 
deviations from orthogonality, indeed |5'pg(0)| < 0.05 for 
all p ^ q, independently of T. Instead, orthogonality is 
clearly violated for strong intramolecular barriers. Off- 
diagonal terms can take values of even a 60% percent 
of the diagonal ones. Moreover, deviations are enhanced 
by decreasing temperature. Fig. |BJ: shows results for the 
unnormalized diagonal terms C'pp(O) (see above). In the 
low p-range the data can be described by an effective 
power law C'pp(O) ~ For fully-flexible chains we 

find approximate gaussian behavior, C'pp(O) ~ [28j . 
However, the introduction of internal barriers leads to 
strong non-gaussian behavior. On increasing the barrier 
strength, the effective exponent x increases up to a value 
of 3.8 for the stiffest case, (A'sji^x) = (35,4), at the 
lowest T. The most local effects of the intramolecular 
barriers are manifested by flattening of C'pp(O) at large p. 

The trends observed for intrachain static correlations 
have their dynamic counterparts. Fig. [7] shows the re- 
laxation times Tp, of the normalized mode correlators 
*I'pp(*) = Cpp{t)/Cpp{Q), as a function of the mode in- 
dex p. We display data for several values of the bend- 
ing and torsion constants {Kb, Kt) and temperatures T. 
The relaxation times have been operationally defined as 
<I>pp(Tp) = 0.3. Data can be again described at low-p by 
an effective power-law Tp ~ The observed trends 

are analogous to those found for the static correlations 
(Fig. [BJ:). Rouse behavior {x = 2) is observed only in the 
fully-flexible limit. Again, as for the static amplitudes 
C'pp(O), X is weakly dependent on T [s^l but strongly 
dependent on the barrier strength, taking higher values 
for stiffer chains. The a;- values for Cpp(O) and Tp at the 
same {K'q,K^) and T are similar. This suggests that 
the structural origin of the observed dynamic anomalies 
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FIG. 8: Normalized mode correlators $pp{t) for stiff chains 
with {Kb,Kt) = (35,4). Panel (a): simulation results at 
T = 1.48. Panel (b): MOT numerical solutions at T = 0.63. 
In both panels, the absolute time is rescaled by the relaxation 
time ri of the p = 1 mode. 



is mainly controlled by intrachain static correlations. 

Fig. shows simulation results for the normalized 
mode correlators for {Kb,Kt) — (35,4), at 

T = 1.48. Times are rescaled by the relaxation time 
of the first mode, ri. Several salient features are re- 
vealed. First, the unambiguous presence of a long-time 
plateau for the modes p = 3 and p = 5, followed by 
an ultimate slow decay. It must be stressed that this fea- 
ture is not related to the structural a-relaxation. Indeed, 
the plateau arises at times far beyond the a-time scale 
(tq, ~ 5 X 10~'^ri for the considered T). This feature is 
instead intimately connected to the relaxation of the in- 
ternal torsional degrees of freedom of the chain. Indeed 
we observe (not shown) that for fixed bending constant 
Kb, the long-time plateau tends to vanish as the value 
of the torsional constant Kt is decreased. 

The observed long-time plateau constitutes a clear 
breakdown of the Rouse model, which predicts single, 
purely exponential decays of the mode correlators (see 
above) . Its origin can be temptatively understood as fol- 
lows. The relaxation of the pth-mode is equivalent to the 
relaxation of a harmonic oscillation of wavelength N/p . 
In the case of strong torsional barriers, the wavelengths 
of some particular modes probe characteristic lengths 
over which chain deformation involves a strong energetic 



penalty (due to the presence of the barriers). Thus, at the 
time scales for which the barrier amplitudes are probed, 
the relaxation of such modes becomes strongly hindered, 
leading to the observed long-time plateau regime and ul- 
timate slow relaxation. Another intriguing feature of 
Fig. [B^, also inconsistent with the Rouse model, is the 
non-monotonous p-dependence of the mode correlators 
at intermediate times prior to the long-time plateau (see 
data for p > 4). 

Now we demonstrate that all the former dynamic 
features can be rationalized in terms of the PRISM- 
based MCT approach of Chong et ai. As exposed in 
Refs. [l^HQij the MCT equations for the unnormalized 
Rouse correlators Cpq{t) are derived as the g — > limit 
of the equations for the self site-site density correlators 
F^Aq,t). The equations for Cpq{t) read 



insitycor 

Mm- 



kBT 



SopSoq + 



kBT 



m 



EpkCkqit) + 



k=0 



kBT 



k=0 



J2 / dt'mpkit~t')Ckqit') = 0, (12) 



with Cpq{0) the inverse matrix of Cpg(O), and Epq = 

{l-So^p>C,)Cpq'{0)/N^,. 

The memory kernel is given by mpq{t) = 
(p/Gtt^) / dkk^S{k)c\k) PrpFlj (k, t)P,qf{k, t) 

[51I - I53I . Thus, prior to solve Eq. (fT2|) . we obtained 
the density-density correlators f{k,t) and self site-site 
density correlators Ffj{k,t) from their respective MCT 

equations (see Ref. [131 )• The static quantities Cpq{0) 
and Cpq^{0), which also enter Eq. (fT2l) as external inputs, 
were directly computed from the simulations at the 
respective lowest investigated temperature. 

Fig.[8]shows a comparison, at er ~ 0.2, of the MCT so- 
lutions for the normalized mode correlators ^pp{t) [panel 
(b)] of the stiffest investigated chains, with the respective 
simulation results previously discussed [panel (a)]. A full 
correspondence between MCT solutions and simulation 
trends is obtained. These include the long-time plateaux 
for p = 3 and p = 5, as well as the sequence in the 
complex, non-monotonous p-dependence for p > 4 at in- 
termediate times. As previously done for the simulation 
data, we can obtain the theoretical relaxation times Tp 
from the condition $pp(Tp) = 0.3 in the theoretical corre- 
lators. The p-dependence of the simulation and theoreti- 
cal times are compared in Fig. [71 at common 0.2, for 
several values of {Kb,Kt)- Again, MCT solutions are in 
semiquantitative agreement with the anomalous trends 
of simulations, with similar exponents for the effective 
power-laws. 

As we observed in Fig. [6] for '^pq{0), there are off- 
diagonal terms of the intrachain static correlations which 
are non-orthogonal. This non-orthogonality persists over 
long time scales, as can be seen in Fig. [SI The lat- 
ter shows simulation and theoretical results for normal- 
ized Rouse cross-correlators ^pq(t), with p = 3 and 
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FIG. 9: Normalized cross-correlators ^pq{t) (for fixed p — 3) 
of stiff chains with {Kb,Kt) = (35,4). Panel (a): simulation 
results at r = 1.48. Panel (b): MOT numerical solutions at 
T — 0.63. In both panels, the absolute time is rescaled by the 
relaxation time ri of the p = 1 mode. 



FIG. 10: Normalized bond correlators Ph{t) for several values 
of the barrier strength at selected temperatures (see legend). 
Panel (a): simulation results. Panel (b): MOT numerical 
solutions. In both panels the absolute time is rescaled by the 
relaxation time n of the p — 1 mode 



q — 1,3,5,7,9. Data correspond to the same tempera- 
tures and barrier strength (the stiffest investigated case) 
of the diagonal correlators of Fig. |8l Again, MCT quali- 
tatively reproduces simulations trends for the case of the 
off-diagonal terms. 

Finally, it is worth noting that the good agreement 
between simulations and MCT for the Rouse correla- 
tors is similar for other observables probing chain dy- 
namics. The reason is that, through the transforma- 
tion Xp(t) = ^^^i Pjp^j{t) (see above), such observ- 
ables can be expressed in terms of the Rouse diagonal 
and cross-correlators . An example is given by the ori- 
entational bond correlator Ph{t) = (b(0) • h{t)) / {b"^ {0)) , 
where h(t) is the bond vector joining two consecutive 
monomers. Following the former transformation we find 

(b(o) . hit)) ^ Ef=r'E;!^=o [^pTi+i - ^p:/][^«>i - 

Pg'j]{Xp{t) ■ X,(0)), where P"^ is the inverse of the 
matrix of coefficients Pjp. Note that this expression 
is exact (the Rouse model makes the approximation 
(Xp(t) • X,(0)) = for p ^ g). Since MCT solutions pro- 



vide the Rouse correlators for all {p, q), insertion of these 
in the former exact expression directly provides Ph{t). 
Fig. [TU] shows simulation and MCT results of Ph{t) for 
several values of {Kb,Kt:) from the fully-flexible limit 
to the stiffest investigated chains. As in previous fig- 
ures, times are rescaled by the respective ri, and data 
in both panels correspond to a common separation pa- 
rameter er ~ 0.2. MCT reproduces semiquantitatively 
the observed simulation trends. These include, on in- 
creasing barrier strength, a relative speed up and slowing 
down (in terms of the seated time t/ri) of respectively the 
short-time and long-time dynamics. MCT also accounts 
for the emergence, for strong barriers, of a plateau at 
t/ri ~ 10~^ and a change in the concavity of the decay. 

IV. DISCUSSION 

On Section III we have shown that, concerning the crit- 
ical temperature Tc, MCT reproduces qualitative simu- 
lation trends for low and moderate barriers. However a 
strong disagreement is found on approaching the limit 
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of stiff chains. We have also found a clear discrepancy 
in the trends of the A-exponent, with a nearly constant 
value from theory and strongly barrier-dependent values 
from simulations. In this Section we discuss possible ori- 
gins of these discrepancies. 

A. Three-point static correlations 

In Ref. [131 we showed that the failure of the MCT 
predictions for strong intramolecular barriers was not 
apparently related with the breakdown of the PRISM 
approximations, which are invoked in the derivation of 
the MCT equations for polymers. Indeed the quality 
of such approximations appeared to be the same for all 
the range of barrier strength here investigated by simula- 
tion and MCT. Despite the mentioned discrepancies be- 
tween theory and simulation, the phenomenological anal- 
ysis of simulation results in terms of a huge set of general 
asymptotic laws of MCT was consistent [13, H^. This 
means that the dynamic exponents involved in the dif- 
ferent tested laws could be, in each case, related to a 
same A'^^. As we discussed in [2^, such scaling laws 
are a mathematical consequence of the bilinear depen- 
dence of the memory kernel on the density correlators 
(see Eq. (O). The specific numerical values of A (and 
by transformation, of the other dynamic exponents) are 
determined by the static quantities entering the vertex 
([H]) 23, 43]. Given the consistency of the phenomenolog- 
ical analysis we speculated that, by retaining the bilinear 
form of the MCT memory kernel, there may be missing 
static contributions in the vertex which are not signifi- 
cant for low barriers, but become increasingly important 
as the limit of stiff chains is approached. Including them 
and solving the MCT equations accordingly, might raise 
the theoretical values of Tc and A, leading to a better 
agreement with the simulation trends. 

Thus, we suggested that intrachain three-point static 
correlations should be explicitly included in the MCT 
vertex. Chain stiffness induces a strong directionality 
in the intrachain static correlations, at least at near- 
neighbor distances. It has been shown that directional- 
ity in static correlations can break the static convolution 
approximation of MCT, Eq. ([6]). A well-known exam- 
ple is given by silica, a network-forming system. For the 
latter the inclusion of three-point static correlations in 
the MCT vertex significantly improves the comparison 
between theory and simulation, with respect to the solu- 
tions obtained under the convolution approximation [sj] ■ 

The calculation of the three-point static correlations 
involved in Eq. (jS]) is very demanding. This is because 
most of the computational time is consumed by the inter- 
chain three-point correlations. For intrachain three-point 
correlations the computation is not demanding. Fortu- 
nately, in the present case only the latter is necessary, 
since the directionality of correlations is only relevant 
along the chain. Thus, the convolution approximation 
is retained for interchain correlations, and it is modified 
only to include the intrachain three-point correlations. 



With this, the new MCT vertex reads '55'! 

V(q, q - k) - ^S{q)S{k)S{\<l - k|)[q • kc(fc) 

+q • (q - k)c(|q - k|) + pq^c^iq, q - k)]^, (13) 

where C3(q, q— k) is the three-point intramolecular direct 
correlation function, given by 



P^C3(q,q-k) = 1- 



^3(q,q- k) 
uj{q)uj{k)u}{\(i~k\)'' 



(14) 



and W3(q, q — k) is the three-point intramolecular struc- 
ture factor 



W3iq,q- 



NrN„ 



E E X 

1=1 a,b,c=l 

exp{i[-q • r^ + k • r^ + (q - k) • r^]}. 



(15) 



Indices in ([T5|) are defined as in ©. The convolution ap- 
proximation for intrachain correlations assumes 0^3 (q, q— 
k) = w(g)a;(fc)aj(|q— k|), or equivalently C3(q, q— k) = 0, 
reducing the vertex (fT3)) to the original Eq. (|8]). 

Figs. [11] and [T^] show representative tests of the con- 
volution approximation for respectively fully-fiexible and 
stiffest investigated chains. Following the scheme pro- 
posed in Ref. [20| , the vectors q, k and p = q — k define 
the sides of a triangle, the first two enclosing an angle 
4> given by cos(j> = (g^ -I- — p^)/2qk. Panels (a) and 
(b) in Fig. [TT] show a test of the corresponding expres- 
sion for an equilateral triangle, uj^{q,q,q) = ui^{q). Pan- 
els (a) and (b) in Fig. [T^] show a similar test for equal 
moduli k ~ q and all the relative orientations (given by 
cos (p) of q and k. In other words, we test the approxima- 
tion uj3{q,q,p — qy/2{l — cos0)) — uj'^{q)u}{p). Data in 
Fig.HHare represented as a function of cos (f) for two char- 
acteristic wave vectors, corresponding to the first mini- 
mum and second maximum of the respective co^^q^q^q) 
(see Fig.HH). 

As already noted in Ref. [13] , the convolution approxi- 
mation for intrachain static correlations provides a good 
description of ujs in the fully-flexible limit. As expected, 
the quality of the approximation decreases by introduc- 
ing intramolecular barriers. Still it constitutes a good 
approximation for all the investigated barrier strength. 
In the case of wave vectors around the first peak of S{q), 
?max ~ 7, the quality is almost unaffected by the bar- 
rier strength, i.e., the terms C3(qmax; qmax — 

k) will be 

small even for the stiffest investigated chains. It must be 
noted that the MCT kernel is usually dominated by the 
contributions around (/max- Thus, the former observa- 
tions suggest that the inclusion of the three-point static 
correlations will modify weakly the MCT solutions ob- 
tained under the convolution approximation. We con- 
firm this by obtaining numerical solutions with the ver- 
tex (|13p . for which we compute the input quantities in- 
volved in Eqs. p3[T5|) directly from the simulations. The 
so-obtained values of the critical temperature Tc and A- 
exponents raise by ~ 1% as much, even for the stiffest 
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FIG. 11: Comparison between tJ3(q,g, q) (dotted lines) and 
the convolution approximation oj''^{q) (solid lines). Panel (a); 
fully-flexible chains at T = 0.50. Panel (b): Stiffest investi- 
gated chains at T = 1.22. 



FIG. 12: Comparison of u)3,{q,q,p — q^2(l — cos 0)) (lines) 
with the convolution approximation Lj'^{q)Lj{p) (symbols), for 
two selected values q — 4.5 and 7.5. Panel (a): fully-flexible 
chains at T = 0.50. Panel (b): Stiffest investigated chains at 
T = 1.22. 



chains, with respect to the previous values (Fig. ^ found 
under the assumption C3 = 0. With all this, we conclude 
that the observed discrepancies between simulation and 
theoretical trends of Tc and A are not related to the break- 
down of the convolution approximation for static three- 
point correlations. The latter indeed retains its validity 
for all the investigated range of barrier strength. 

B. Dynamic heterogeneities 

It is well-known that the quality of the Kawasaki ap- 
proximation for dynamic correlations (see above) breaks 
on decreasing temperature. This feature is specially crit- 
ical around the a-time scale [36l - l40j , leading to the com- 
plete failure of the MCT predictions associated to it, 
as the power law behavior D~^,Ta ^ (T — Tc)^^ , or 
the time-temperature superposition of density correla- 
tors. This breakdown is usually assigned to the emer- 
gence of strong dynamic heterogeneities on approaching 
the glass transition [SGl - lioj . Having noted this we may 



speculate that, for some reason to be understood, in- 
creasing the barrier strength strongly enhances dynamic 
heterogeneities. This might result in a lower quality of 
the MCT and might be the reason for the observed dis- 
crepancies between theory and simulation trends for Tc- 

Now we show that this is not actually the case, and 
that there is no correlation between barrier strength and 
enhanced dynamic heterogeneity. Non-gaussian param- 
eters provide a simple way of quantifying the strength 
of the dynamic heterogeneity. They display large posi- 
tive values at the time scales for which the respective van 
Hove function strongly deviates from the gaussian limit. 
This occurs when a significant fraction of particles has 
performed displacements very different from the average. 
Here we discuss the two most popular non-gaussian pa- 
rameters. The standard or 'fast' parameter is given by 
a2{t) = f((Ar(t))4)/((Ar(t))2)2 - l. The 'slow' param- 
eter introduced by Flenner and Szamel [56| is defined as 
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flexible (circles) and stiff chains with {Kb,Kt) — (25,1) 
(squares). Empty and filled symbols correspond respectively 
to the fast {a2{t)) and slow (72(i)) parameters. 



72(t)-i((Ar(i))^)(l/(Ar(t))2)-l. ^ 

By construction Q2 (t) and 72 {t) are identically zero for 
a gaussian form of the van Hove self-correlation function. 
However, as noted in [sl], large positive values of these 
parameters have a very different microscopic origin, re- 
flecting distinct aspects of dynamic heterogeneity. In the 
case of the fast parameter a2{t), large values originate 
from a significant population of particles which have per- 
formed much larger displacements than the average. This 
effect is generally maximum at the time scale t* around 
the end of the caging regime. Thus, 02 (t) increases from 
zero at i = up to a maximum at t*, and decays to zero 
at longer times. The increase of the maximum a2{t*) 
on decreasing temperature reflects a progressive enhance- 
ment of dynamic heterogeneity, at the decaging process, 
on approaching the glass transition. 

The slow parameter 72 (i) exhibits analogous trends 
for the temperature and time-dependence. However the 
maximum of 72 (t) takes place at much longer scales than 
i*, namely around the a-relaxation time Tq. This effect 
originates from a significant population of particles which 
at t ~ Tq have performed much smaller displacements 
than the average [s^. 

Fig. [T^ compares simulation results of 02 (i) and 72(^)1 
for the fully-flexible case and for very stiff chains. For a 
fair comparison we have selected temperatures at which 
the respective a-relaxation times are similar. These are 
T — 0.50 and T = 1.05, for respectively fully-flexible 
and stiff chains, and correspond to a separation param- 
eter CT ~ 0.04 (see Fig. [3|). Let us remind that the a- 
time scale can be estimated, e.g., as /(gmax,0 = 0.2. 
In Fig. [T3]we display /(gmax,^) for both systems at the 
former temperatures, showing that the respective a-time 
scales are roughly the same, Tq ~ 10"*. The decaging 
times, which can be estimated from the start of the de- 
cay from the plateau in /(ftnaxiO: ^^^^ roughly the 
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FIG. 14: Density-density correlator /(gmax, t) for fully-flexible 
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same, t* ^ 500. This equivalence is indeed reflected in 
the trends of the non-gaussian parameters in Fig. 1131 
Thus, in both systems a2{t) is peaked at t* ~ 500 and 
72 (t) is peaked at Tq ~ 10"*. 

Having noted this equivalence of time scales, data in 
Fig. [13] do not reflect any enhancement of dynamic het- 
erogeneity on increasing the barrier strength. Actually, 
the opposite effect is suggested by the lower values of 
a2(t) and 72 (i) for stiff chains with respect to the fully- 
flexible case. With this, we discard a major role of dy- 
namic heterogeneities as the reason for the observed dis- 
crepancies between simulations and MCT solutions for 
very stiff chains. 

C. Chain packing 

As mentioned in Section III and shown in Fig. [5l MCT 
fails at reproducing the peak around qc ^ ^ for the q- 
dependence of the KWW times of density-density corre- 
lators. As noted in Ref. [2Cy for the fully-flexible case, 
the origin of this peak may be related to dynamic corre- 
lations between centers-of-mass of the chains. The latter 
might arise from the effective packing between the poly- 
mer coils, interacting as fully penetrable spheres of size 
i?g. This interpretation is not clear in view of the results 
of Fig. O since the value of qc does not seem to be related 
with 27r/i?g. Having noted this, Chong et al. found that 
the incorporation of the static correlations between the 
centers-of-mass in the MCT equations did not improve 
the description of the simulation results. As shown in 
[20t . this is not unexpected due to the almost featureless 
form of the static structure factor of the centers-of-mass 
S'cm('7)- The inset of Fig. [5^ shows simulation results of 
ScM{q), for the same barrier strength and temperatures 
of the KWW times in the main panel. The introduc- 
tion of chain stiffness does not induce significant features 
in ScM{q), appart from a stronger signal at low q. The 
latter indeed suggests that packing effects between the 
polymer coils are even weaker than for the fully-flexible 
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case. Within the former interpretation, this would be 
consistent with the lower intensity of the mentioned peak 
of /Tq^^^ at qc ~ 4. All these results suggest that dis- 
crepancies between simulations and MCT on increasing 
chain stiffness are not related to a dynamic coupling, not 
accounted for within the theory, to the slow modes at 
qc ~ 4. Indeed, this coupling seems to be weaker for stiff 
chains. 

D. Outlook 

In summary, in this section we have discussed possible 
origins for the discrepancies, concerning the structural 
relaxation, between simulations and MCT on increasing 
barrier strength. We discard a major role, in comparison 
with the fully-flexible case, of three-point static correla- 
tions, dynamic heterogeneities and chain packing. These 
effects become even weaker on increasing chain stiffness. 
We remind that such effects are indeed neglected in the 
derivation of the MCT equations used here (see Section 
III). Results in this section suggest that this is not less 
justified for very stiff chains than for fully-flexible ones. 

How to improve the theory to account for dynamic 
trends in stiff chains is an open question. A way might 
be the reformulation or extension of the MCT equa- 
tions, retaining the bilinear form of the kernel, in terms 
of new dynamic observables coupled to density fluctua- 
tions. Such observables can be adequate for describing 
particular dynamic features which are not captured by 
the usual observables, i.e., the number density fluctua- 
tions p[(\,t). Some examples are roto-translational site 
fluctuations adapted to the molecular symmetry, as has 
been shown, e.g., for dumbbell-like molecules [iTl, |57| or 
for a simple model of orthoterphenyl [l^ . The inclusion 
of density fluctuations of centers-of-mass improve results 
for rigid molecules |l8l| concerning a peak in T^/Tq^^^ 
at intermediate q |58| , similar to that observed here at 
qc ^ 4. As discussed above, this is not the case for poly- 
mer chains. Though there is no characteristic symmetry 
in polymer chains, roto-translational density fluctuations 
can also be defined over sites a, b at some characteristic 
distance |a — 6|, perhaps probing the relevant length scale 
27r/gci which according to the data of Fig. [S] seems to de- 
pend weakly on the barrier strength. Whether this pro- 
cedure may improve the agreement between MCT and 
simulations remains to be solved. 

V.CONCLUSIONS 



By means of simulations and solution of the equations 
of the Mode Coupling Theory, we have studied the role of 
intramolecular barriers, of arbitrary strength, on several 
aspects of polymer dynamics. The investigated dynamic 
range extends from the caging regime characteristic of 
glass-formers to the relaxation of the chain Rouse modes. 
Solutions of the MCT for the structural relaxation re- 
produce qualitative trends of simulations for weak and 
moderate barriers. However, a progressive discrepancy 
between MCT and simulations is revealed as the limit 
of stiff chains is approached. We have tested the validity 
of several assumptions inherent to the theory. Deviations 
from the theoretical predictions do not seem related with 
dynamic heterogenities, which indeed arc not enhanced 
by increasing the barrier strength. Moreover, the convo- 
lution approximation for three-point static correlations 
retains its validity for stiff chains. Even the role of slow 
modes at intermediate length scales, not accounted for 
by MCT, becomes less significant on increasing chain 
stiffness. At this point it is not clear how to improve 
the MCT equations in order to remove the mentioned 
discrepancies for the case of stiff chains. We have sug- 
gested the possibility of formulating the MCT equations 
in terms of roto-translational density fluctuations over 
specific length scales. 

Concerning the relaxation of the chain degrees of free- 
dom, MCT provides a microscopic basis for the observed 
deviations from the Rouse model on increasing the bar- 
rier strength. These include anomalous scaling of re- 
laxation times, long-time plateaux, and non-monotonous 
wavelength dependence of the mode correlators. Be- 
yond usual phenomenological models for chain dynamics 
(the Rouse model being the corresponding one for fully- 
flexible chains), MCT provides a unified microscopic pic- 
ture down to time scales around and before the a-process, 
which is not accounted for within the mentioned models. 
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